Quantum critical end point of the Kondo volume collapse 
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The Kondo volume collapse describes valence transitions in /-electron metals, and is characterized 
by a line of first order transitions in the pressure-temperature phase plane terminated at critical end 
points. We analyze the quantum critical end point, when the lower end point is tuned to T = 0, 
and determine the specific heat, thermal expansion, and compressibility. We find that the inclusion 
of quantum critical fluctuations leads to a novel bifurcation of the first order phase line. Finally, 
we show that critical strain fluctuations can cause both, superconductivity and non-Fermi liquid 
behavior near the critical point. 
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Quantum criticality in heavy fermion materials is pre- 
dominantly discussed in the context of magnetic phase 
transitions 0. In contrast, instabilities in the charge sec- 
tor, which are well known to occur in Ce and Yb based 
intermetallics, are usually first order transitions. In the 
case of strong first order transitions, like the Kondo vol- 
ume collapse (KVC) transition between a and 7 Ce 
or the related behavior in YbInCu 4 [3} .critical fluctua- 
tions are irrelevant. However, theories [J, |5j and exper- 
iments for the KVC yield a line of first order phase 
transitions in the pressure-temperature phase plane that 
is terminated at critical points on the high and (some- 
times) on the low temperature end. The location of this 
lower end point is tunable. For example, a suppression of 
the volume collapse transition was achieved by alloying 
Ce with Th or byapplying an external magnetic field 
to Ceo.sLao.iTho.i [TJ. In the case where the lower end 
point is tuned to T = 0, a quantum critical end point 
emerges. 

In this paper we investigate the behavior in the vicinity 
of the quantum critical end point of the Kondo volume 
collapse. As shown in Fig. 1, we demonstrate that the 
p— T phase transition line is qualitatively changed by crit- 
ical fluctuations. In addition, we determine the pressure 
and temperature dependence of the heat capacity, the 
compressibility, and the thermal expansion, and discuss 
the possibility of superconductivity and non- Fermi liquid 
behavior caused by these critical fluctuations. Ce and Yb 
based intermetallic alloys with a moderate high temper- 
ature bulk modulus are good candidates to identify such 
a quantum critical point of the Kondo volume collapse, 
where the critical fluctuations are in a universality class 
similar to that of a metamagnetic end point 0|. Being 
in the charge sector, KVC fluctuations may serve as a 
mechanism for s-wave superconductivity 0, . Strong 
retardation effects of the critical pairing interaction are 
important and enhance the superconducting transition 
temperature as compared to a BCS calculation. Our 
predictions for the temperature dependence of the com- 
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FIG. 1: Pressure-temperature phase diagram for the KVC 
model for different values of the bare bulk modulus B from 
the slave boson theory (see text). The main figure includes 
critical fluctuations, the inset is the mean field result. Lines 
represent first order transitions from a low Tk phase, to the 
left of the line, to a high Tk phase, to the right of the line 
(this is for Ce alloys, the reverse would be true for Yb ones). 
The solid circles and dot indicate critical end points for T > 
and T = 0, respectively. The inclusion of fluctuations causes 
the first order line to "bifurcate" for B < B* . The lower line 
then terminates at a quantum critical end point for B — B* . 



pressibility and the thermal expansion are unique for the 
quantum critical end point of the Kondo volume collapse. 
Identifying such a behavior in a superconducting mate- 
rial strongly suggests that the pairing is due to critical 
strain fluctuations. Finally, recent advances in the the- 
ory of quantum phase transitions in metallic svstems|ll| 
demonstrate that previous results obtained for magnetic 
critical points are unstable against non-analytic correc- 
tions to the collective mode dynamics. Quantum criti- 
cality in the charge sector, as discussed in our paper, is 
therefore one of the very few cases where robust predic- 
tions for non-Fermi liquid behavior can be made. 

A candidate for this behavior is CeCu2Si2- a ;Ge a ; under 



2 



pressure A volume discontinuity at p ~ 15GPa was 
observed in CeCu2Ge2 13']. Evidence for a valence in- 
stability in CeCu2Si2 at pressure p ~ 25kbar was already 
given in Ref.Q The increase of the residual resistivity 
and the drop of the T 2 coefficient of the low temperature 
resistivity with pressure in CeCu2Si2 and CeCu2Ge2 have 
been interpreted in terms of critical valence fluctuations 
0,0]. Also, the ability to separate two superconduct- 
ing regimes by varying pressure in CeCu2Sii.2Geo.s [13 
supports the point of view that two different mechanisms 
are at work, one related to the magnetic critical point and 
the other to the valence fluctuation one |lCj. 

The order parameter at a volume collapse transition 
between two isotructural states is the trace 



tre = 



V-Vn 



To 



(1) 



of the strain tensor e, where Vo is a reference volume. 
Our choice for Vo as the volume at the quantum critical 
end point will be discussed below. The mean field theory 
for the KVC transition 0, [j| starts from the Gibbs free 
energy 

G( P ,T) = F K (e,T)+p^ + W (2) 



where B is the bare (high temperature) bulk modulus 
that describes the elastic properties in the absence of the 
Kondo effect at the volume Vo. In a cubic system B = 
\ ( c ?i + 2cJ 2 ) with bare elastic constants and c\ 2 , P is 
the pressure, and Fk (e, T) the Kondo contribution to the 

free energy. For T « T K , F K (T, V) ~ -V (t k + ag) 

with Kondo temperature Tk — De~ x / JK (D is of order 
the bandwidth of the conduction electrons) and a of or- 
der unity. Due to the exponential variation of Tk with 
Jk , moderate variations of the Kondo coupling Jk with 
volume cause a strong volume dependence of the Kondo 
temperature. This, in turn, yields a non-linear equation 

of state p (V) = — §|f| t = — Pq 7%?\t- ^ ^ s eas y to ver ify 
that with appropriate values of B, even a linear volume 
dependence of Jk yields p — V isotherms similar to those 
of the van der Waals theory of the liquid-gas transition 
;4t, 5] , with a first order transition where a discontinuous 
volume change occurs. For sufficiently large B, the line 
of first order transitions is terminated at a lower critical 
end point. For T 3> Tk, Fk is dominated by a saturat- 
ing spin entropy, resulting in a termination of the line of 
first order transitions at a high T critical end point |l6| . 

In order to go beyond mean field theory and to analyze 
the quantum critical fluctuations of the KVC transition, 
we consider e (x,r) as a space and time dependent fluc- 
tuating strain with a field theory governed by the action 



S [e] = Sa [e] + d d x dr ($ (e (x, r)) + pe (x, r)) . 



Here, Sa [e] describes strain fluctuations that are nonlocal 
in time and space while the static potential $ (e) includes 
the nonlinear effects of the theory. The external pressure 
p occurs in the last term such that the Gibbs free energy 
follows as G (p) = — TTog J deexp(-S). We first derive 
Eq. |3 from a microscopic model. Then we discuss the 
physical consequences that follow from the analysis of 
S[e}. 

To describe the Kondo lattice physics of the system, 
we start from the infinite U Anderson lattice model 



(4) 



where Hq = CkC^Cko- is the Hamiltonian of the con- 
duction electrons, e° the bare / electron energy, and t 
the hybridization. The model describes strongly corre- 
lated electrons because of the infinite Coulomb repulsion 
of the / electrons, leading to n fi = J2a=i fLf™ < 1 ( N 
is the degeneracy of the / electrons). 

In a situation where the value of the hybridization de- 
pends sensitively on the volume of the system, i.e. on the 
local strain e i; we replace H — » H + H c with 



(7 f 



2 



(5) 



where 7 is the coefficient of the assumed linear volume 
dependence of the hybridization. B is, as in Eq. |5J the 
bare bulk modulus and Nq the number of unit cells in 
the system. In Eq. we assumed that e only couples to 
the hybridization. Additional couplings to the conduc- 
tion band and / level energies also occur, but will not 
qualitatively change the conclusions of our paper. They 
will be ignored in what follows. 

In order to obtain a quantitative insight into the low 
temperature behavior of the Kondo lattice, we use a slave 
boson mean field approach [TtI ] to enforce the constraint 
n fi < I- The main results of this paper do not depend 
on the details of this method, though. The result of the 
slave boson mean field theory is the free energy $ (e) per 
volume Vo as a function of strain. Within mean field 
theory $ (e) = Fk/Vq + \Be 2 of Eq. and results for 
the volume collapse transition are in full agreement with 
Refs.iandli 

In Fig. 2 we show the results for the pressure-strain 
dependence of the Gibbs free energy G(e) — <&(e) + pe 
for different values of the bare bulk modulus as obtained 
from the slave boson mean field calculation. Close to the 
transition and at low T the potential has the form 



<f> (e) = $0 -p*s 



3 



(6) 



(3) 



Here, $0 is the e independent part of <5(e) and p* (T = 0) 
the pressure at the critical point. The coefficients p* , b, 
v and u can be obtained from the slave boson theory by 
differentiating the free energy with respect to the strain. 
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FIG. 2: Slave boson mean field results for the strain depen- 
dence of the Gibbs free energy. A double minima indicates 
the presence of a first order transition. For B — B* and 
T = 0, the double minima collapse to a single one, indicating 
the presence of a quantum critical end point. For B < B* 
(inset), the transition is always first order. 



They are temperature dependent due to the temperature 
dependence of the spin entropy, and thus at low tempera- 
tures are functions of ^- , where T is of order the Kondo 
temperature. Note that b also has a term determined by 
the bare bulk modulus B (Eq. 2). The lower end point of 
the KVC transition is located at T = 0, if ff = f^f = 0. 
This is achieved by varying B to the value B* and the 
pressure to p* . In order to have this transition take place 
at e = 0, we choose the so far arbitrary reference vol- 
ume Vo such that <I> (e) has no contribution of order e 3 
at T = 0. This implies that v = xT 2 /T with yt > and 
that p* now depends on B. 

The dynamic and nonlocal part S& of the action can 
also be derived within the slave boson approach. Inte- 
grating out the fermionic degrees of freedom, 



(n(g)-n(o))| £ (g)r 



fq 



where J q ... = TJ2 n J 

n ( q ) = -NjH 2 J q , ((;■;;:(;!/ , s + cri c;;, : „ 



with q 



(7) 

(q,<x>„) and 

, The polar- 
ization bubble II (q, to) determines the dynamics of the 
strain field and is determined by the single particle propa- 
gators of the fermions within the slave boson calculation, 
i.e. G^ c (r) = - (T r c q (r) c\ (0)), etc. The evaluation of 
this expression for small q and ui (but with vf |q| ~S> to) 
yields^ 




\2k F J 



(8) 



As long as the system is in a heavy fermi liquid regime, 
the generic ui and q-dependence of Sn in Eq. [S] is inde- 
pendent of the details of the slave boson approach. The 



latter does however allow one to determine the prefactors 
F and a in terms of microscopic quantities^! 0|. The 
model, Eqs. and is essentially the same as the one 
used by Millis et al. jg] in the context of metamagnetic 
quantum criticality in metals and numerous analogies ex- 
ist. 

From the effective action, Eqs.|21Eland|Sl we obtain the 
equation of state including fluctuation corrections and 
determine the location of the volume collapse transition. 
It follows that 



p = p* - (b+ u(e 2 ))e + 



ve 



ue 



(9) 



with (e 2 ) 



qq 

2k F 



There are 



' B=B" ~ Jq \ v F \q\ 

two distinct sources of temperature variation in Eq. [5] 
First, due to the temperature dependence of the mean 
field entropy, it holds for T < T : p* (T) = p* + 

b(T) = b — Cj^, etc., where £, C are positive con- 
stants. Second, the order parameter fluctuations are 
temperature dependent: (e 2 ) T = (£ 2 ) T=0 + C,T^~ with 
i2.i5r 1/3 



c 



> for d — 3. For d < 5, the latter 



a S / 3 (-U J? fc F ) 1 /3 

is the dominant effect (fluctuation corrections due to the 
cubic term cx v are subleading). 

We first ignore the T dependence of (e 2 ), formally cor- 
responding to d > 5, and analyze the mean field p — T 
phase diagram. The critical end point is located at T = 
if p = p* and B - B* = b + u(e 2 ) T=Q = 0. Close 

to this transition at T = 0, e (B = B*,p) oc (p* - p) 1/3 
and e{B,p = pl) oc (B* - B) 1 / 2 . At finite T we ob- 
tain a first order transition if B < B* + t^rr- For 

' J-a 

B > B* + C^- ; no transition takes place. The zero tem- 
perature transition is first order for B < B* , second order 
for B = B*, and absent for B > B* . For the following 
set of parameters: t = 0.5eV, D — leV, = — 2.34eV, 
(yielding T K ^ 0.012D and n f = 0.87), B* = 254kbar, 
Pq = — 4.7kbar and 7 = 2.14, the corresponding p — T 
curves are shown in the inset of Fig. 1. The value for 7 is 
chosen to yield results that quantitatively reproduce the 
mean field theory of Ref.Q. 

The phase diagram changes qualitatively if we include 
the temperature dependence of the critical fluctuations 
of the order parameter (e 2 ^ T , which is dominant at low 
temperatures for d < 5. Besides entering with a different 
power than the mean field terms, the key new aspect of 
this fluctuation effect is that it leads to an increase, as 
opposed to a decrease, of the b (T) + u (e 2 ) term in Eq.EI 
as T increases. While the mean field entropy softens the 
bulk modulus, critical fluctuations harden it. Thus, if 
the critical point is located at T = 0, increasing temper- 
ature will suppress the transition rather than making it 
go first order. This effect is demonstrated in Fig. 1, where 
we show the fluctuation corrected transition lines for dif- 
ferent bulk moduli. The competition between the fluctu- 
ation term (which dominates at the lowest temperatures) 



4 



and the mean field terms (which dominate at higher tem- 
perature) leads to the development of a gap in the first or- 
der line at a B > 0.96835*. With increasing B, the first 
order line below this gap eventually collapses to a quan- 
tum critical point at B — B* , and then disappears for 
larger B. Thus, the critical end point becomes isolated in 
the p — T phase plane. For the parameters we used, the 
gap in the first order transition ~ 0.002Z? ~ 20K should 
be observable. 

Next, we use the equation of state Eq. El close to the 
quantum critical end point to determine the pressure and 
temperature dependence of e, the compressibility k and 
the thermal expansion (3. The analysis of Eq. [5] yields 



follows that 



oc 



Ih 



AV 
Vdf 



]_dV 
V~dp 



ocT" 



(10) 



where 



1 



and 8 = | for d > 3 and 



7 



and 6 = ^ 



We obtain the exponents 

3 ' ' 

for d < 3. For d < 3 the T-dependence of the strain is 
governed by the fluctuation corrections in the equation 
of state, while for d > 3, the mean field T dependence 
■yr- in <f> (e) of Eq. [fj] dominates. For d — 3 both terms 
yield exactly the same T dependence. The strain fluc- 
tuations also give rise to a singular heat capacity Hf- 
oc -logT for d = 3 and ^ oc T^r for d < 3. For 
the physically relevant heat capacity at constant pressure 



C P = C V + VkT (dp/dT) z v we find C p 



C v oc T 3 -t, i.e. 
the dominant contribution to C p comes from Cv- Pres- 
sure tuning at T — yields a diverging compressibility: 
k oc \p — Pol 2 ^ 3 - Critical fluctuations in the charge sec- 
tor furthermore yield non-Fermi liquid behavior of the 

3-d 

single particle self energy: T,k F (w) oc uj \uj\ 3 for d < 3 
and Sfc F (a;) oc Lulog \u>\ for d = 3. 

Finally we comment on the possibility of superconduc- 
tivity caused by the interaction of electrons with these 
critical fluctuations. Strain fluctuations couple to elec- 
trons just like phonons or other bosonic charge excita- 
tions, leading to an attractive interaction in the s-wave 
pairing channel. Superconductivity as caused by valence 
fluctuations was discussed in Refs. © and OH In those 
treatments, retardation effects of the pairing interaction 
enter the theory solely via the upper cut off, wq, of the 
theory. In the weak coupling limit this leads to the 
well known BCS-formula for the transition temperature 



rpBCS 



uj e 



where A„ 



is the dimensional 



coupling constant of the pairing interaction. However, in 
the vicinity of the critical end point, the strain fluctua- 
tions serve as pairing boson which becomes massless. As 
shown in Ref. Il9t I2C1 retardation effects are more subtle 
in the case of pairing due to a massless boson governed 
by Eq. |SJ For a given value of X p , the transition tem- 
perature is enhanced compared to the BCS result and it 



T, 



uj e 



(11) 



Thus, the critical fluctuations discussed in this paper are 
good candidates for bosons that cause s-wave supercon- 
ductivity. 

In summary, we have analyzed the critical behavior in 
the vicinity of a T = end point for the Kondo volume 
collapse transition. We find diverging specific heat and 
thermal expansion coefficients, and a diverging compress- 
ibility. We also find that critical fluctuations break apart 
the line of first order transitions, leading to a novel "isola- 
tion" of the critical end point in the pressure-temperature 
phase plane. Finally, we discussed the implications of our 
results for the existence of superconductivity near a KVC 
quantum critical end point. 
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